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ABSTRACT 

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization 
of molecular layers above the continuum forming photosphere. 

Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the 
spectral domain 1.45-2.50 fim with a spectral resolution of « 35 and baselines ranging from 15 to 88 meters. We perform independent 
image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model. We also compare the 
data to ID dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic 
giant branch (AGB) stars. 

Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots 
whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At 
« 2.00 /um and in the region 2.35 - 2.50 //m, the photosphere appears extended and the radius is larger than in the H band. In this 
spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular 
layers. The wavelength dependence of the visibility can be qualitatively explained by ID dynamical models of Mira atmospheres. The 
best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of = 8.82 + 0.50 
mas. The H 2 molecule seems to be the dominant absorber in the molecular layers. 

Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmso- 
pheres. In particular, we see molecular (water) layers that are typical for Mira stars. 

Key words, stars: AGB, Post- AGB, supergiant - stars: atmospheres - stars: individual: VX Sgr - techniques: interferometric 



1. Introduction 

VX Sagittarii (HD 165674) is a coo l semi-regular variabl e 
with a long mean period of 732 days (iKholopov et all Il987b . 
lLockwood & Wingl (1 1982b reported a spectral type varying from 
M5.5 (near the time of visual maxim um) to M9.8 (at mini- 
mum light). lLockwood & Wingl (Il982|) determined that the ef- 
fective temperature of VX Sgr is ran ging between 3300 and 
2400 K (maximum to minimum light). iGarcia-Hernandez et aD 
( 2007) find r e ff =2900 K at the time of their high-resolution spec- 
troscopic observation, when the star was near minimum light 
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(AAVSCQ). lLockwood & Wingl d 1982b found that VX Sgr ex- 
hibits stronger CN and VO bands with respect to Mira variables 
with similar temperature. Enhanced CN absorption is an indi- 
cator of high luminosity i n RSGs of e a rlier t ype and, together 
with VO, also of S stars. ISpeck et alJ (12000b categorized VX 
Sgr as an oxygen-rich star and found a strong silicate feature at 
10 yum, that indicates a dusty circumstellar environment. Using 
aperture-m asking and IR/optical- telescope array interferometry 
at 2.16 //m jMonnier et al.1 (12004b revealed that VX Sgr exhibits 
a dusty environment with a flux contribution of about 20% in the 
K band and some evidence of departure from circular symmetry, 
even if they could not place strong limits on possible asymme- 
tries because of calibration unc ertainties. The dusty e nvironment 
is confirmed by HST images dSchuster et all l2006h . VX Sgr's 
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circumstellar environment is the result of the heavy mass loss 
experienced by the star (1.3 • 10~ 5 Moyr _1 , CO measurements; 
KnaprjetaLL Il989l) . The mass-loss proces s appears to be par- 
ticula rly asymmetric for the inner regions (IChapman & Cohenl 
11986b - Using AAVSO data, iKamohara et al.l (120051) show that 
the optical light curve has a much smaller amplitude of about 
2 mag in the years 1998-2003, much less than the usual 6-7 
mag. An examination of AAVSO data shows that the decrease 
to this smaller amplitude has happened several times in the last 
70 years, and that the star is probably currently in that state. 

The classification of VX Sgr as a red supergiant or an AGB is 
thus not firmly established. A further constraint can be brought 
by an estimate of its luminos ity, in order to b e tter a scertain 
its position in the HR diagram. iHumphreys et al.l (Il972|) placed 
VX Sgr in the vicinit y of the Sgr OBI cluster at 1.7 kpc; 
iMurakawa etafl (120031) fou nd 1.8+0.5 kp c measuring the GHz 
H2O maser expansion; IChen et' al. (2007) reported a distance of 
1.57+0.27 kpc using 43 GHz SiO maser pr oper motions; finally , 
the trigonometric parallax of Hipparcos (Ivan Leeuwenl 120071) 
gives a distance of 0.262(+0.655/-0.109) kpc, probably unreli- 
able due to the size and asymmetry of the stellar photosphere. 

Using AAVSO data, we find VX Sgr was at maximum lu- 
minosity during our observations, and we as sume a T e s of 320 
to 3400K. With the 2MASS K magnitude dCutri et all l2003h . 
assuming a dista nce d= 1 .7 kpc, and usin g data for Galactic red 
supergiants from iLevesque et al.l (120051) . we derive a luminos- 
ity logL/L = 5.25 + 0.25 (Af bo i = -8.4 + 0.6). The error 
bar accounts for uncertainties in the photometry, and in the as- 
sumed r e ff at the time of our observation, impacting the bolo- 
metric correction at K. Circumstellar emission in IR may in- 
crease further the luminosity by a few tenths of magnitude. 
Putting the star at 1.3 kpc would decrease the estimated lumi- 
nosity to logL/L© = 5.00 + 0.25. The radius is then about 
1200 R Q , and A v =2 to 4. This is a too hig h luminosity for 
an AGB star (e.g.. IVassiliadis & WoodL fl993|) . Even compared 
to so-called super-AGB stars, where the most recent models 
show a maximum of log L/L©«_4. 8 with typical masses rang- 
ing between » 7 — 11 M© (ISiessL l2006h and logL/L© « 5 
dPoelarends et alll2008l) , VX S gr's luminosity is extremely high. 
iGarcia-H ernande z et al.l (120091) found AGB stars with similar lu- 
minosities and masses of ^6-7 M© showing Rb enhancement 
in the Magellanic Clouds, and they argue that these AGB stars 
may be more lumi nous due to a contributio n from Hot Bottom 
Burning. However, IGarcia-Hernandez et al.l (120071) . using syn- 
thetic spectra based on classical hydrostatic model atmospheres 
for cool stars with extensive line lists, found VX Sgr to be the 
only not Li-rich, long-period, high OH expansion velocity star 
of their Galactic AGB sample. On the other hand VX Sgr's low 
effective temperature, and large V variabil i ty are quite untyp- 
ical for an RSG, although ILevesque et al.l (120071) found high- 
variab ility, low-r e ff RSGs in the Magellanic Clouds. Illeger et al] 
(1 19971) studied the pulsations properties in red supergiants from 
10 to 20 M Q with high luminosity to mass ratio and show that 
very large pulsation periods, amplitudes and mass-loss rates may 
be expected to occur at and beyond central helium ex-haustion 
over the time- scale of the last few 10 4 years. This could lead to 
an overall dimming of the star after a period of stronger oscil- 
lations subsequent to enhanced mass-loss and ejection of a dust 
shell that screens the stellar radiation. 

It appears that the evolutionary status of VX Sgr is still not well 
established, and more investigations are needed. In particular its 
chemical composition should be scrutinized. 

We discuss here interferometric observations of VX Sgr 
made with the VLTI/AMBER instrument in the near IR. The aim 



of this paper is to study the continuum forming photosphere, the 
existence, and the characterization of molecular layers of VX 
Sgr probing different wavelengths in the H and K bands. 



2. Observations and data reduction 

We obtained near-infrared interferometry data of VX Sgr 
with the Very Large Telescope Interferometer (VLTI, 



Haguenauer et al, 2008) using the near-infrared beam combiner 
AMBER dPetrov et all 120071) that simultaneously covers the 
J, H, and K bands with a spectral resolution of « 35. VX 
Sgr has been observed in less than 1.5 months using the AT 
configurations: A0-D0-H 0, D0-H0-G1 an d E0-G 0-H0. The 
fringe tracker FINITO dLe Bouquin et all 12008b has been 
used for all the observations. In addition to the science target, 
three calibrator stars have been observed close in time and 
interleaved with VX Sgr: HD169916 (K0IV), HD146545 
(K5III) and HD166 295 (K2III/IV) The calibrator diameters 
were retrieved from IRichichi et al.l (120051) . iBorde et al.l (l2002h 
and IRichichi & Percheronl (l2002k . The diameter errors are 
of the order of 1%. Details are reported in Table [T] with the 
projected baseline lengths (B^), the position angles (PA^), the 
spectral interval, the calibrators used, the optical seeing, and the 
coherence time. 

Raw visibilities and closure phases were comp uted with the 
latest version of the amdlib package (release 2.2, ITatuUi et all 
120071) and the yorick interface provided by the Jean-Marie 
Mariotti Center. Individual frames were averaged after select- 
ing the 20% best frames based on the fringe SNR only and with 
a piston smaller than ±20 /im. We decided to discard J band 
fringes because the data quality was significantly worse than for 
longer wavelengths. 

In addition we used the addons by iMillour et al.l (120081 re- 
lease 1.53) to calibrate the datasets. The transfer function mea- 
surement for one night is shown in Fig.[T] Both science and cal- 
ibrator targets have the same detector integration time of 0.05 
ms, and after the calibration, VX Sgr data have been averaged. 
The error bars on the calibrated visibilities include the statistical 
error of averaging the single frames, the errors of the calibra- 
tion stars' angular diameters, and the scatter of the transfer func- 
tion measurements. This scatter (top panel of Fig. [B is much 
larger than the internal errors, computed by the reduction soft- 
ware. This basically means that the visibilities errors (between 
0.05 and 0.1) we use in this paper reflect calibration issues af- 
fecting simultaneously a whole range of wavelengths, while the 
wavelength-to-wavelength error is much smaller (typically 0.01 
to 0.05): while the visibility errors seem large in a single dataset, 
the wavelength- variation errors of the same dataset are small. 

The absolute wavelength correction has been done using 
the telluric Kitt Peak spectra which we convolved to match the 
spectral resolution of the AMBER data. In VX Sgr data, the 
band gaps (i.e., between J and H and H and K) are visible and 
we made a linear two-component adjustment of the wavelength 
scale which gave a systematic offset of -0.21 jim with respect to 
the initial AMBER table and a 7% wavelength stretch. 

Fig. [2] shows the final UV-plane coverage of all observa- 
tions that successfully passed all steps of the data reduction and 
calibration quality control. The north-west south-east direction 
is not completely covered because of the actual AT geometry, 
while the east-west direction is particularly favored. The star 
exhibits large wavelength-dependent visibilities and has a clear 
non-zero, non-180 closure phase (bottom panel of Fig. []]), evi- 
dencing asymmetries in the intensity distribution. 
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Table 1. Observations log for the AMBER observations of VX Sgr. All observations were carried out using FINITO and an integra- 
tion time of 50 ms. 



Date 


Bp 


PA P 


Spectral 


Calibrators 


Seeing 


Coherence 




[m] 


[°] 


range \jim] 




["] 


time [ms] 


2008-05-24 


H0-D0 (58.76)/D0-A0(29.37)/A0-H0 (88.14) 


-99.8/-99.8/-99.8 


1.8731-2.4766 


HD169916 


0.65 


4.5 


2008-05-26 


G1-D0 (71.55)/D0-H0(63.94)/H0-G1 (71.48) 


-46.0/70.5/7.0 


1.4492-2.3398 


HD 166295 


0.97 


6.0 


2008-06-06 


H0-G0 (30.10)/G0-E0(15.06)/E0-H0 (45.17) 


-119.0/-119.0/-119.0 


1.4492-2.4458 


HD 1699 16, 


0.79 


4.5 










HD 146545 






2008-06-07 


H0-G0 (28.64)/G0-E0(14.33)/E0-H0 (42.97) 


-98.8/-98.8/-98.8 


1.4655-2.5000 


HD169916 


1.16 


2.0 


2008-06-08 


H0-G0 (31.98)/G0-E0(16.00)/E0-H0 (47.99) 


-108.0/-108.0/-108.0 


1.4454-2.3708 


HD 1699 16, 


0.60 


2.5 










HD146545 






2008-07-03 


H0-G0 (29.63)/G0-E0(14.82)/E0-H0 (44.46) 


-100.0/-100.0/-100.0 


1.4578-2.2786 


HD 1699 16, 


1.29 


2.0 










HD146545 






2008-07-04 


G1-D0 (63.13)/D0-H0(35.39)/H0-G1 (68.78) 


-4.3/35.4/68.8 


1.4454-2.4062 


HD169916 


0.41 


7.5 


2008-07-05 


G1-D0 (71.52)/D0-H0(63.92)/H0-G1 (71.46) 


-46.1/63.9/6.9 


1.4511-2.4832 


HD169916 


0.65 


6.5 


2008-07-06 


G1-D0 (65.26)/D0-H0(55.29)/H0-Gl (71.48) 


-55.2/52.6/-7.78 


1.4540-2.5211 


HD169916 


0.60 


3.5 
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Fig. 2. UV-coverage for of the AMBER observations of VX Sgr. 
The radial extension of the uv tracks reflects the spectral cov- 
erage of our AMBER measurements, covering the H and the K 
bands. 



Time (h) 

Fig. 1. Transfer function for the night 2008-07-03. Top panel: 
squared visibilities averaged over the region 2.1-2.2/mi for both 
the calibration stars (grey with error bars) and the science star 
(green with error bars). The computed transfer function is ap- 
proximated by a linear fit (solid black line), and dashed lines 
display the scatter. Bottom panel: closure phase averaged over 
the same spectral region for both the science target (red) and 
calibrators (black). 



3. Image reconstruction 

The first step in our analysis is a chromatic image reconstruc- 
tion of our data to probe different layers in the photosphere and 
above. 

The image recons truction process is similar to the one per- 
formed for T Lep in iLe Bouquin et al.l (120091) . For the recon- 



struction, we use th e MIRA software package (iThiebauti [2008: 
ICotton et all 120081: ILe Besnerais et all l2008h . The image is 
sought by minimizing a so-called cost function which is the sum 
of a regularization term plus data-related terms. The data terms 
enforce agreement of the model image with the measured data 
(visibilities). The regularization term is a x 2 minimization be- 
tween the reconstructed image and an expected image. The ex- 
pected image is issued from a preliminary image reconstruction, 
strongly constrained by the assumption of circularity. Each spec- 
tral bin has been processed independently (Fig.O left column). 
The reconstructed images clearly highlight different behaviors 
across the wavelength range: in the H band (« 1.45 - 1.80 yum), 
the intensity distribution is inhomogeneous and a bifurcation of 
the image core into a few bright "spots" is visible; at « 2.00 
/mi and at the upper K band edge (« 2.35-2.50 jim) the radius 
appears extended and much larger than in the H band. Artifacts 
may be introduced by the poor UV-plane coverage in one direc- 
tion (north-west south-east, see Fig [2]). However, the detection 




Fig. 3. Left column: reconstructed images of VX Sgr for several AMBER spectral bins across the H and K bands. The intensity, /, 
is normalized to the range [0, 1] and plotted as I 033 . Some contour lines are indicated (0.251, 0.275, 0.314, 0.372, 0.588, 0.889, 
0.929, 0.976,0.983,0.990,0.997). The resolution of the interferometer is illustrated in the bottom left part of each image by the PSF 
of a 88x70 m telescope. Right column: images representing our best-fitting geometrical model for the same spectral channels. The 
intensity, /, is normalized to the range [0, 1] and plotted as I 033 . 
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of inhomogeneities is out of doubt because there are clear signa- 
tures in the closure phases (Fig.Q]). 

4. Geometrical toy models 

The reconstructed images are not good enough to firmly estab- 
lish the presence of inhomogeneities on VX Sgr's photosphere 
because the different sources of errors can smooth out the infor- 
mation. Thus, in order to confirm their presence and to constrain 
their flux relative to the total flux, we used a model-fitting ap- 
proach. We approximated the star with a uniform disk and its 
extended shells by a Gaussian disk, and we added a series of 
point sources to model the spots. We performed a fit using all the 
observed wavelengths together, but used five regularly spaced 
reference wavelengths to effectively compute the model (1.54, 
1.78, 2.03, 2.28, and 2.52 /mi). To match the data with the model, 
we interpolated the parameters for other wavelengths using cu- 
bic spline interpolation. The global optimum of the fit was found 
using a set of simulated annealing algorith ms complemented by 
standard gradient descent algorithms (see iMillour et all 120091 
for a first application of this approach). A model with three spots 
(i.e., performing a fit with a uniform disk, a Gaussian disk, and 
three spots) significantly enhance the fit (x 2 ~ 25), especially 
for the visibilities and closure phases up to spatial frequency of 
« 240 cycles/arcsec. Spatial frequencies above this value corre- 
spond to a single measurement at the largest baseline length. The 
details about the fits are reported in Table [2] 

Our toy model is probably not perfectly describing the object 
at the highest angular resolution available, but gives a resonable 
fit to intermediate angular resolutions. We also note that two of 
the bright spots are located at the position of the stellar disk, and 
one has to be outside the stellar disk in order to fit the observed 
closure phases. We estimate that the upper value of the spots' 
flux is < 10% of the total flux and that the flux of the spot lo- 
cated outside the photosphere has very likely a flux equal to zero 
in the H band. Fig. [3] (right column) displays the appearance of 
the toy model at some selected wavelengths. As a final remark, 
these images resemble the reconstructed images in the left col- 
umn of the figure. We also note that the bright spot located out- 
side the photosphere is particularly visible in the lower right part 
of each panel of the reconstructed images at lower wavelength 
and somehow less evident at longer wavelength even if it is ex- 
pected by the geometrical model. 

What we can conclude from this model-fitting approach is 
the following: 

- a model with three spots can explain the visibilities and 
closure phases for spatial frequencies < 240 cycles/arcsec. 
Higher spatial frequencies are probed in our dataset by only 
one observation and indicate a more complicated object 
shape at smaller scales; 

- The best-fitting model locates two of the spots inside the sup- 
posed photosphere of the star, but the third spot, mandatory 
to explain the closure phase deviations, is located outside the 
photosphere and it is brighter at longer wavelenths; 

- The spots have individual fluxes which, at all wavelengths, 
do not exceed 10% of the total flux of the star, for all wave- 
lengths. 

5. Model fitting 

The visibility data of VX Sgr show significant wavelength- 
dependent features. We fitted the visibility curves and the closure 
phases using a Gaussian disk for each spectral channel, which 




H 2 Q H 2 Q H 2 Q - 



1.4 1.6 1.8 2.0 2.2 2.4 
Wavelength (|am) 

Fig. 4. VX Sgr Gaussian FWHM values as a function of the 
wavelength compared to the model Ml 8 predictions. The dashed 
line indicates the small est values. The position o f the H 2 bands 
are also indicated (after lLancon & WoodLl2000l) . 



gives more robust results than a uniform disk model. Figure 01 
illustrates the change of apparent size with wavelength. This 
is consistent with what is visible in the reconstructed images 
(Fig. [3]): the radius gets larger between 1.8 and 2.1 yum, reach- 
ing local maximum at 2.0 yum, then gets smaller again to a mini- 
mum around 2.15 to 2.25 /mi, eventually it becomes much larger 
toward 2.5 fim. However, Fig. 01 gives only a rough idea of the 
apparent diameter of VX Sgr because the intensity profile is not 
a Gaussian disk. 

To reproduce the wavelength d ependence, w e used th e one- 
dimensional dynamic models of llreland etall (2004a : b) for 
oxygen-rich Mira stars that include the effect of molecular layers 
in the outer atmosphere. These are self-excited dynamic mod- 
els whose grey atmospheric temperature stratification was re- 
computed on the basis of non-grey extinction coefficients that 
contain all relevant molecular absorbers (e.g., H 2 0, CO, TiO; 
solar abundances) but do not contain dust. The stellar parame- 
ters were assumed to be close to those of the M-type Mira proto- 
types o Cet and R Leo. This model ser ies has been successfully 
used bv lWittkowskietal.1 (120071 l2008h to explain AMBER ob- 
servation of the Mira star S Ori. Due to the poorly understand- 
ing of the VX Sgr stellar parameters, there is no fully consistent 
model of this star, and the models we used can show only typical 
characteristics of such a star. Figure \5\ shows the comparison for 
short (left and central panels) and long (right panel) projected 
baselines. The wavelength dependence of the visibility is similar 
for all the nights. 

For the fit, we fol lowed the fitting procedure used by 
IWittkowski et al.l (120081) and used only the short baseline data 
obtained with E0-G0-H0, because we do not expect to match 
details of the intensity profile probed at high spatial frequencies 
with the M series that has not been designed to match a star 
like VX Sgr . The dynamic models predict large changes in the 
monochromatic radius caused by the molecular layers above the 
continuum-forming region in the stellar photosphere. The band- 
passes around 1.9 and 2.5 /mi are significantly affected by H2O 
molecules with some contribution of CO in the H band and in 
the long-wavelength part K band. The mole cular-band effects 
chang e significantly with phase as shown in IWittkowski et al.l 
(120081) . The model Ml 8 (see Table[3]) provides the best fit to the 
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Table 2. Parameters of the best geometrical toy model (x 2 ~ 25). 



A 


Gaussian disk 


Uniform disk 




Point source 




Point source 






Point source 






size 


flux 


size 


flux 


X 


y 


flux 


x y 


flux 


X 


y 


flux 


[jim] 


[mas] 


[%] 


[mas] 


[%] 


[mas] 


[mas] 


[%] 


[mas] [mas] 


[%] 


[mas] 


[mas] 


[%] 


1.54 


<20 


26.7 


8.3 + 2.5 


67.0 


0.6 + 0. 


8 2.8 + 0.5 


<10 


-0.7 + 0.7 -0.5 + 0.5 


<10 


10 + 4 


-8.5 + 3.5 





1.78 


< 20 


32.1 


8.1 + 1.5 


59.3 






<10 




<10 









2.03 


<20 


47.3 


8.0 + 0.6 


43.8 






<10 




<10 






<5 


2.28 


< 30 


46.1 


7.7 + 0.4 


47.5 






<10 




<10 






<10 


2.52 


28.7 + 26 


66.9 


12.4 + 1.9 


22.9 






<10 




<10 






<10 





Fig. 5. Measured VX Sgr visibility data (red crosses with error bars) compared to the visibilities derived from: (i) the ID Mira 
atmosphere model Ml 8 (model with a phase O = 0.75), black solid line; (ii) the best-matching synthetic image of the 3D simulation 
of a RSG star (top central and right panel in Fig. [6]), light blue triangles; (iii) best match synthetic images of the 3D simulation of 
an AGB star (bottom row in Fig. [6]); blue dots. Left and central panels belong to the night 2008-06-07, and the right panel is for the 
night 2008-07-03. 



VX Sgr data among the 20 available phase and cycle combina- 
tions of the M series. 

The M model qualitatively explains the overall wavelength 
dependence of the apparent stellar radius (Fig. solid black 
line). However, while the fit is very good at short baselines, the 
comparison is less good at long baselines. Two reasons can ex- 
plain this discrepancy: (i) the stellar parameters of the M models 
are not well suited for VX Sgr, (ii) there are some surface inho- 
mogeneities detected in the data, that are not in the M model. 
The angular photospheric diameter corresponding to the refer- 
ence radius at 1.04 jim is ©mis = 8.82 ± 0.43 mas for short base- 
lines where the fit is better. The error on the diameter includes 
systematic calibration and model uncert ainties. This result is i n 
agreement to what has been found by Monnier et al. I d2Q04h . 



©2.16/zm = 8.7 + 0.4 mas . However, our models do not include 
dust while iMonnier et all found a dusty environment around the 
star with a flux contribution of about 20% in the K band. 

We also fitted Gaussian FWHM values to the synthetic visi- 
bility values based on the Ml 8 model intensity profiles for each 
spectral channel, using the same procedure that was used to fit 
the measured data. The comparison between the fit to the Ml 8 
model intensity profiles and the fit to the measured data is shown 
in Fig. ffl 

The available Mira models show good agreement with VX 
Sgr's data. However, this comparison can just give a basic qual- 
itative picture of VX Sgr's surface and more detailed interpreta- 
tions of these data must be addressed to next generation models 
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Table 3. Parameters of the simulations. All the models have solar 
metallicity. 



Model 


Stellar 


M 


L 


R 




logg 


name 


type 


[M Q ] 


[L Q ] 


[R Q ] 


[K] 


[cgs] 


M18^ 


Mira 


1.2 


4840 


210 


3310 




st28gm06n06^ 


AGB 


1 


6935 


429 


2542 


-0.83 


st35gm03n073 


RSG 


12 


93 000 


832 


3490 


-0.34 



a ID mode l; the radius i s the R osseland radius R=0.81R P at Mira 

phase ^0.75; (Ireland et al., 2004ah 

b 3D model; Freytag and Hofner (2008) 
c 3D model; Chiavassa et al. ( 20091) 



with stellar parameters appropriate for VX Sgr, which are cur- 
rently not available. 

6. Complementary comparison with 
three-dimensional simulations 

In this section, we quantify if the observed asymmetries are con- 
sistent with three-dimensional simulations of surface convec- 
tion in RSG and AGB stars, and if their chaotic photospheric 
structure is adequate to explain the different surface layers struc- 
tures observed The simulation are carried out w ith C0 5 BOLD 
dFrevtag et all 120021: iFrevtag and Hofnerl 120081) . Parameters of 
the models are given in Table [3j The pulsations are not artifi- 
cially added to the models (e.g. by a piston) but are self-excited. 
Excitation by not- stationary sonic convective motions are re- 
sponsible for the pulsations. Molecular opacities are taken into 
account, but radiation transport is treated in grey approximation, 
ignoring radiation pressure and dust opacities. Dynamical pres- 
sure lets the density drop much slower than expected for a hy- 
drostatic atmosphere. In our 3D simulations, the average density 
drops exponentially, and there is no sign of a wind or an ex- 
tended shell with relatively lar ge densities. Some more tec hni- 
cal information can be found i n l Freytag and Hofnerl (120081) . the 
CO 5 BOLD Online User ManuaQ, and in a forthcoming paper by 
Freytag (2010, in prep.). 

We computed « 3500 synthetic images from these simula- 
tions at the same wavelengths of the observations using our 3D 
radiative transfer code OPTIM3D dChiavassa et alll2009h . Then, 
we generated v isibility curves using the method described by 
IChiavassa et all In the case of the RSG simulation, the synthetic 
images do not change strongly in diameter, shape and asymme- 
tries across the wavelength range of the observations, and, hence, 
the observed wavelength dependence of the visibility is not ex- 
plained by this approach (blue triangles in Fig. [5]). The synthetic 
images of the AGB simulation show a noticeable wavelength 
dependence and their agreement with the observed visibilities is 
slightly better than the RSG model, in particular at short base- 
lines (left and central panel, blue dots of Fig. [5]). Nevertheless, 
the two 3D simulations analyzed do not provide a better fit to 
our AMBER data than the Mira models described in Sect. [5j 

At last, in order to explain the detected spots on the surface 
of the reconstructed images (Fig. [3]), we compare the synthetic 
images corresponding to the best-fitting visibility curves of the 
3D simulations (see Fig. [5]) to the reconstructed image at 1.6 yum 
(Fig [6]). This wavelength corresponds roughly to the H" 1 con- 
tinuous opacity minimum (i.e., the photosphere becomes more 
transparent). The reconstructed image shows two spots on the 

2 www. astro . uu . se/~bf/co5bold_main . html 




-4-2 2 4 
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Fig. 6. Reconstructed image at 1 .6 fim with the intensity normal- 
ized to the range [0, 1] and plotted as I 033 . The contour lines are 
the same as in Fig. [3] 

surface. Figure [7] displays the synthetic images of the RSG star 
which show large convective cells 4 - 5 mas), and on top 
of that there are small-scale granules (« 0.5 - 1 mas). The ob- 
served visibility points do not show a good agreement especially 
at long baselines where the expected signal is much lower than 
the simulation's predictions. 

The AGB star (Fig. [7]) displays a pulsating stellar atmosphere 
with strong asymmetric structures caused by the shocks; the con- 
vective envelope is hardly visible because of the optical thick- 
ness of the at mosphere even without the inclusion of the dust 
opacity tables dFreytag and Hofneil 120081) . The AGB's visibil- 
ities and closure phases show large departures from symmetry 
already in the first lobe, but even in this case, the long baseline 
data points cannot be explained. 

The presence of spots on the reconstructed image of VX Sgr 
can be qualitatively explained by the synthetic images, even if 
the distribution of the spots in our simulations has little chance 
to be exactly the one we observed with AMBER. Unfortunately, 
we cannot constrain the 3D simulations in term of surface in- 
tensity contrast using these observations because we could not 
determine an accurate value for the spots' flux level (Table [2]). 

7. Conclusions 

Our AMBER observations unveil, for the first time, the shape of 
VX Sgr's surface. 

The individual wavelength image reconstruction has been 
carried out using the MIRA software. VX Sgr displays visibili- 
ties and images with strong wavelength dependence: (i) surface 
asymmetries in the H band (« 1.45 - 1.80 yum) and (ii) an ex- 
tended radius at « 2 and 2.35-2.50 yum. We find that a geometri- 
cal toy model composed of a uniform disk plus a Gaussian disk, 
and three spots gives a reasonable fit for the AMBER data. We 
claim that at least two spots are present on the photosphere of 
VX Sgr, and we show that this is qualitatively predicted by 3D 
hydrodynamical simulations of stellar convection. In addition, 
the toy model shows that one spot resides outside the photo- 
sphere and it is brighter at longer wavelengths. The presence 
of this spot located outside the photosphere could be related to 
a complex and irregul ar structure in the surr ounding of VX Sgr 
as already detected bv lKervella et al.l (120091) for the RSG alpha 
Ori. Also the two bright spots which appear at the position on 
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Fig. 7. Top left panel: synthetic image from the RSG simulation at 1 .6 yum, convolved with a 5.9x4.6 mas PSF. The intensity is plotted 
as in Fig.O The simulation has been scaled in size to match approximately the visibility data points. Top right panel: original RSG 
simulated image, the range is [0; 300 000] erg cm" 2 s" 1 A" 1 . Bottom panels: synthetic visibilities and closure phases derived from 
the above image (black solid line and dots) compared to the observations (red with error bars) at 1.6 yum. The dashed lines indicate 
the minimum and maximum amplitude of variations of the visibilities and closure phases issued from different rotations of the 
image. In bottom right panel, the axisymmetric case is represented by the solid lines. 



the photosphere may both originate at the depth of the contin- 
uum photosphere, or higher up in the molecular layers. In fact, 
these spots appear at all wavelengths (Fig. [3]), also at the water 
bands, which hide a large part of the continuum photosphere. 
This might be a hint that they originate in the molecular layers 
far above the continuum photosphere. 

We used ID dynamical oxygen-rich Mira model predictions 
to explain the visibility data points, and we found that the fit is 
very good at short baselines. We conclude that H2O molecules 
strongly affect the visibility, and thus this molecule seems to be 
a dominant absorber in the molecular layers. The atmospheric 
structure of VX Sgr seems to qualitatively resemble Mira-star 
models which show molecular layers above the continuum form- 
ing layer. In addition, its photosphere shows bright spots that 
can be related to giant-cell surface convection and possibly to 
the molecular layer. However, we must point out that the Mira 
model used here has stellar parameters that are not consistent 
with what we expect for VX Sgr, and there are only observa- 
tions from a single epoch to do the interpretation. Due to various 
uncertainties, the classification of VX Sgr must be further in- 
vestigated in the near future by, firstly, combining multi-epochs 
photometric and spectroscopic observations to better determine 
its stellar parameters; and secondly, by observing with interfer- 
ometers at higher spectral resolution to study in greater detail the 



changes across the H and K bands using new theoretical models 
with consistent stellar parameters. 
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